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STATE-VECTOR CONTROL APPLIED TO LATERAL STABILITY 


OF HIGH-PERFORMANCE AIRCRAFT 

By Waldo I. Oehman and Jerrold H. Suddath 
Langley Research Center 


SLM4AKY 




State-vector control has been applied to the problem of lateral stability 
augmentation of high-performance aircraft. The problem considered was that of 
evaluating the feedback gains that would produce desired stability character- 
istics. The feedback gains are obtained as solutions of linear algebraic equa- 
tions. A numerical example illustrates the practicability of the method. / y 

ikO 


INTRODUCTION 




Perturbed lateral motions of a high-performance aircraft may have undesir- 
able stability characteristics at some flight conditions. A linear feedback 
control system usually is designed to augment the lateral stability. Ifce design 
problem is quite complex if deflections of either the ailerons or the rudder 
produce changes in both the rolling and yawing moments. Hie cross-control 
effects (yawing moment caused by aileron deflection and rolling moment caused 
by rudder deflection), if large, may require that the feedback of the rolling 
and yawing velocities be applied to the ailerons and rudder to provide addi- 
tional damping. The designer must determine the required cross-damper feedback 
gains (rudder deflection due to rolling velocity and aileron deflection due to 
yawing velocity) that provide desired stability characteristics, in addition to 
the usual yaw-damper and roll -damper feedbacks . 

A well-known mathematical theorem on linear dynamical systems (ref. l) will 
be applied to derive analytical expressions for the gains that will produce any 
desired roots of the characteristic equation of the augmented system. The feed- 
back gains are obtained as solutions of linear algebraic equations if the linear 
dynamical system has the property of being completely controllable vising a 
single linear control function. Hie property of complete controllability has a 
precise mathematical definition, and it is possible to check this property of 
the system and to know in advance whether or not the required gains can be com- 
puted. A numerical example is included for illustration. 


SYMBOLS 




span, ft 


rolling -moment coefficient. 


Rolling moment 
qSb 


g 

*X 

J Z 

j - 


yawing-moment coefficient, 

qSb 

side-force coefficient, — — - - orce 

qS 

constant vector with components c-^ and Cg 

acceleration due to gravity, 52.2 ft/sec^ 

moment of inertia about principal body X-axis, slug-ft^ 

moment of inertia about principal body Z-axis, slug-ft^ 


feedback gain vector with components k^ where i is 1, 2, 
mass, slugs 

dynamic pressure, ipV^, lb/sq ft 

wing area, sq ft 
time, sec 

control vector with components 6 a and 6 r 
airspeed, ft /sec 

state vector with components x-^, X 2 , x^, and x^ 

angle of attack of principal X-axis 
angle of sideslip 
aileron deflection 


3, and k 


5 . 


rudder deflection 



parameter 


* 

>1 roots of characteristic equation where i is 1, 2, 3, an d 4 

p air density, slugs /cu ft 

0 angle of roll 

if angle of yaw 



dp 



= dCj 
^a 





_ dCn 
" ^ 

_ ^n 

d6 a 


n - ^ 

" <** 

°Y ‘T 1 

All angles are in radians unless otherwise noted. One dot over a quantity 
mean s the first derivative with respect to time, and two dots over a quantity 
means the second derivative with respect to time. 


STATE-VECTOR EQUATIONS OF MOTION 


A system of linear differential equations with constant coefficients that 
describe the perturbed lateral motion of an airplane (see ref. 2) may be 
written, using principal body axes, as 
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p - S - 1 * - + ♦ - 0 


$ - SS*. Cjjj - ^ c z 6 a - c z 6 r = °V 

y I Y *p p I Y l6 D a I Y ‘8- r / 


0 


X °a X °r 




c b 


qSb 


Terms involving the rotary derivatives and the effect of control authorities 
on side force have "been omitted for convenience. These terms could have been 
included without affecting the method of analysis, except that computations 
would be more complicated. 

If = 0, x 2 = x 3 = x 4 = then the state of the system 

may be expressed as a four -component vector x^ or 



(: 


Similarly, the control deflections may be written as a two-component vector 
u*, or 


u* = 



( 


The system (eq. (l)), with equations .(2) and (3).> may be written in vector- 
matrix notation, as 



•where 




fi a -1 
V 


0 

0 10 

00- 

I x p 

0 0 0 


<lSb 

0 0 0 


0 

0 


0 

0 

00- 

qsb 

X x z 5 a 

<lSb 

hT c \ 


qSb r 

T 

_z a 

iSb _ 


The elements of the matrices JjTj and [G~j are constant parameters 

on flight conditions and airplane characteristics. The matrix [cf] 
trol effectiveness matrix. 


that depend 
is the con- 


The system of equations (eq. (4)) is the state-vector form of the equations 
that describe the perturbed lateral motion of an airplane with controls. Aug- 
mentation of the stability of the airplane is an a l yzed in the following 
discussion. 


STABILITY AUGMENTATION 


Consider the airplane equations of motion with u* = 0. The system of 
equations of motion for the free system is 

£= [A]5? (5) 

The characteristic equation for the free system is a polynomial equation in A 
and is defined by 
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p(A) = det|Y[Y] - [a] = 0 


or 


p(x) % * 3 + qSb (^ % ' § ■ -») 


0l .ix 2 -fi^cx-o 


( 6 ) 


(7) 


where [Y] is the identity matrix, A is a parameter, and det^^ means 
determinant of a square matrix. The roots Aj_ of the characteristic equation 
provide information about the stability of the airplane. 


If the free airplane system has 'undesirable stability characteristics, -chen 
a feedback control system may be desirable to augment the stability. The prob- 
lem, therefore, is to choose the control vector if such that the augmented 
system (eq. (4)) has the desired stability characteristics, or, equivalently, 
that the roots of the characteristic equation of the augmented system have the 
desired values. 


A mathematic ally simple form of the control vector results if the compo- 
nents 6 a and 6 r are assumed to depend only on a single linear feedback 
function and to vary proportionally. The control vector then can be written 


or 


u = c(k . 3) = c*k* x 


S„ = 


— c i(^i x i ^■ 2 x 2 k 3 x 3 k 4 x 4) 
8 r = c 2 ^k 1 x 1 + k 2 x 2 + + k^x^ j 


(8a) 


(8b) 


In eq-uation (8a), the linear feedback function (k* . x) is the dot product of the 
vectors k* and 5 ?, and k* is the transpose of k*. The vector ? is called 
the feedback gain vector and its components k^ are the feedback gains . The 
components of the constant vector c* are c-^ and C 2 « It is clear from equa- 
tions (8b) that 5 r and 5 a will generally vary proportionately (when c-^ / 0 
and c 2 / 0). However, for augmentation with either rudder or ailerons alone, 
either or may be identically zero with the other component equal to 

one. In any case, it is required that all four state-vector components (p, 0, 

0, and ijf) be fed back. 

The system of equations of motion (eq. (4)), with the control vector 
defined in equation (8a), is 

x = QQx* + [gJc* ^ x* = Qa] + [Yjcf k^jx* (9) 
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Note that 


£»>*] 


]>] + DO 5 *^ 


is the 


is an n x n matrix and that 
matrix of the augmented system. The characteristic equation for tfie augmented 
system (eq. (9)) is also a polynomial equation in A and is 


q(A) = det A[l] - QQ “ 00 = 0 


( 10 ) 


Reference 1 shows that for each choice of characteristic roots, a gain 
vector k can be determined if a condition of complete controllability, as 
defined in reference 3, is satisfied. Furthermore, the gain vector is com- 
puted by solving a set of linear algebraic equations in the gain factors. A 
brief discussion of complete controllability of an augmented system is pre- 
sented in the appendix. 

The determinant of equation (10), when expanded, gives q( A) as 
q (A) = A 4 - (Lckj + N c k 4 + Y p )A 5 

- [oLp - Np + (oLc - N c )ki + Lck2 - YpLckj - YpN c kj^J A 2 


where 


- [f L p + & L c k x - YpL c k 2 + (NpL c - LpN c )k 5 + a(N p Lc - LpN c )k^A 


- (NpL c - LpN c )k 2 - f(N p L c - LpN c )k 4 = 0 



(11a) 


(Hb) 
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The polynomial q(A) also may be written in terms of the desired characteristic 
roots Ax as 

q(?0 = (A - Ax) (A - A 2 )(A - Aj)(A - H) 

= A^ - (Ax + A2 + A^ + A^ + |axA 2 + AjAij. + (Ax ^2) (^5 ^ 


- |AxA 2 (A 5 + A4) + A jAl|_ ( Ax + A 2 jJ A + A-^ApA^A^ - 0 


( 12 ) 


If the coefficients of powers of A in equation (lla) are equated to the coef- 
ficients of like powers of A in equation (12), four simultaneous linear alge- 
braic equations in the gains are obtained. They are 


L c k^ + N c k^ - Ajl + - Y^ 


- (aL c - N c )k 1 - L c k2 + YpL c k^ + Y^N c k^ - + ^ 3 ^% + (^1 + ^2) (^3 + ^ 1 +) + ^p - 

V L c k l ■ Y p L c k 2 + ( N P L C ■ L p N c) k 3 + a ( N p L c ■ L p N c)H = ^ 1 ^ 2 (^5 + \) + \ + Ag) " y L p 

(L p N c - N p L c )k 2 + |(L p N c - W p L c )k^ = 


} (13) 


Equations ( 15 ) are the desired linear equations for the feedback gains kx 
in terms of the parameters of the airplane and the roots Ax of the character- 
istic equation. The ratio Cx/c 2 (when Cx/c 2 = S a /5 r ), with either Cx 
or c 2 equal to 1 , must be specified before these equations can be solved. 
However, it is clear that this ratio should be chosen equal to the ratio of 
augmenter deflection limits, which the designer must set from practical consid- 
erations. Then, the ratio of the control deflections will have the proper 
value when the augmenter system is saturated by a large disturbance. 


EXAMPLE CALCULATIONS 


Consider an airplane having the characteristics and flight conditions 
given in the following table: 
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22.36 
32.2 
5021 
67 199 
390. k 
200 
200 
6000 
125 000 
6 

20 

0.015 

0.31 

- 1.0 

-0.075 

-0.15 

0.08 

-0.108 


The roots of the unaugmented characteristic equation are 

7 ^ = 0 

>2 = 0.00kk9Ck 

= -0.0107835 + 1. 78698 j 
\ = -O.OIO7835 - 1. 78698 j 


b, ft 

g, ft/sec 2 . . 
I X , slug-ft 2 . 
I z , slug-ft 2 . 
m, slugs . . . 

q., lb/sq. ft * . 

S, sq. ft • * • 

V, ft /sec • • • 

Altitude, ft . 
Mach number . . 
a, deg . . * . 

C 1 f per radian 

Cn Q , per radian 
P 

Cy , per radian 

P 

C 1 , per radian 

°a 

Ci , per radian 
b r 

per radian 
Cn^, per radian 


These characteristics are representative of an airplane similar to the X-15 
flying at Mach 6 and an altitude of 125 000 feet. The components of the vec 
tor c*, are 



The selection of values for c-^ and c ^ is based on the assumed relative 

authorities of the ailerons and rudder for augmentation purposes . Maximum 
authorities of the rudder and ailerons should be attained simultaneously for 
the present choice of c*. The desired stability is represented by the roots 
of the controlled system which were chosen to satisfy military handling - 
qualities specifications (ref. 4) and are 
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\ = -0.0546 
7*2 = -0.695 
A 5 = -0.546 - 5-l4j 
\ = -0.546 + 5.l4j 


where j = \J^1. 

The condition for complete controllability, from equation (A7) in the 
appendix, is 


det 



0 


The inequality is satisfied with the numerical data of the example. 

The solution of equation (15), with use of the example data, resulted in 
the following gains: 


kp = 0.959796 
k 2 = 0.0026195 
k 5 = 0.0715657 

k^ = 0.0591162 

A similar set of gains can be computed for any set of desired roots of the 

characteristic equation (eq. (9))> and, in a complete analysis, the roots may 
have to be chosen selectively in order to insure reasonable gains. 

Figure 1 presents the motions of the controlled airplane, for the partic- 
ular set of chosen roots, following an initial sideslip angle of 10. 9° • The 
period of the oscillatory motion is 2 seconds and the time to damp to one-half 
amplitude is 2 seconds. The motions of the free airplane system are included 
in figure 1 for comparison. The free airplane has a period of about 5-5 sec- 
onds, and the time to damp to one-half amplitude is about 64 seconds. Although 
the motions were computed for a large initial sideslip angle, the aileron and 
rudder deflections are within probable limits for an airplane of the type 
considered. 
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CONCLUDING REMARKS 


State-vector control theory has teen applied to a stability problem of 
high-performance aircraft. Analytical expressions have been derived for feed- 
back gains that will produce any desired roots of the characteristic equation 
of the augmented system. The gains are given in terms of the parameters of the 
system and the desired roots. A numerical example illustrates the practica- 
bility of the method. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April J, 1965* 
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appendix 


COMPLETE CONTROLLABILITY 


A linear dynamical system is said to Be completely controllable if, at any 
initial time t, any initial state x can be taken to the origin in a finite 
length of time by the application of a suitable control function. 

Consider the mathematical model of a linear dynamical system 

If. [a]? + [G] u> (Ai) 

where x* is an n- vector, the state of the system, and u* is an m- vector, the 
control of the system. The rectangular matrices |X] and Qj] are, in gen- 
eral, functions of time. However, for present purposes, they are assumed to be 
constant. 

A necessary and sufficient condition for the constant system (eq. (Al)) 
to be completely controllable is 


rank 


00- DGH- D0 2 H- • • • - H” 0 


= n 


(A2) 


Suppose that the control vector u is a scalar s times a known vector 


Define a new n-vector 

w *s [g]X 

so that the dynamical system is represented by 



5? + w*s 


The condition for complete controllability of equation (A5) is 

rank[Mj = rank^, [aJw*, [a] 2 ^, • • • > |X] n ” 1 w*J = n 
or, since QiJ is an n X n matrix. 


(A3) 

(AM 


(A5) 


(A6) 
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APPENDIX 


det 


0*, H 2 ^ 



4 o 


A more detailed and general discussion of 
presented in reference 3* 


complete 


controllability is 
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Figure I.- Concluded. 
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